Ultrafast and spatially resolved studies of charge carriers in atomically-thin 

molybdenum disulfide 
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Atomically-thin molybdenum disulfide is emerging as a new nanomaterial with potential appli- 
cations in the fields of electronic and photonics. Charge carrier dynamics plays an essential role 
in determining its electronic and optical properties. We report spatially and temporally resolved 
pump-probe studies of charge carriers in atomically-thin molybdenum disulfide samples fabricated 
by mechanical exfoliation. Carriers are injected by interband absorption of a 390-nm pump pulse 
and detected by measuring differential reflection of a time-delayed and spatially-scanned probe pulse 
that is tuned to an exciton transition. Several parameters on charge carrier dynamics are deduced, 
including carrier lifetime, diffusion coefficient, diffusion length, and mobility. 
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Molybdenum disulfide is a transition metal dichalco- 
genide with an indirect bandgap of 1.29 eV. 1 It is com- 
posed of monolayers of S-Mo-S that are bound by the 
weak van der Waals force, while the atoms in each layer 
are bound strongly by ionic-covalent interactions. The 
layered structure allows fabrication of atomically-thin 
films where the quantum confinement can significantly 
modify the electronic and optical properties. Although 
few-layer— ~— and even monolayer— M0S2 samples have 
been fabricated and studied since 1960's, they have be- 
come topics of great interest again, partially owing to the 
progress in studies of graphene^ 

Very recently, strong photoluminescence (PL) was ob- 
served from few-layer M0S2 samples*^ which was at- 
tributed to an indirect-to-direct band gap transition 
that occurs when varying the thickness from bulk to 
monolayer*^ Such a transition was also confirmed by 
theoretical calculations 1 ^ and scanning photoelectron mi- 



probe technique to study charge carriers in few-layer 
M0S2 samples fabricated by mechanical exfoliation. Car- 
riers are injected by direct interband absorption of a 390- 
nm pump pulse, and detected by measuring the differen- 
tial reflection of a 660-nm probe pulse. We found that 
the excitonic absorption is reduced by the carriers, which 
is consistent with the phase-state filling effect. By spa- 
tially and temporally resolving the signal, we deduce a 
carrier lifetime of 100±10 ps and a carrier diffusion coeffi- 
cient of 20±10 cm 2 /s, corresponding to a mobility of 800 
cm 2 /Vs and a diffusion length of 450 nm. The demon- 
strated technique can be used as a powerful tool to study 
charge carriers in various MoS2-based structures. 

Few-layer M0S2 samples are fabricated by mechanical 
exfoliation^ with an adhesive tape from natural crystals 
(SPI Supplies). By using a silicon substrate with a 90- 
nm Si02 layer, flakes with different atomic layers can 
be readily identified by using the contrast of microscopy 



croscopy measurements^ In addition to the possible images^^ Figure [T^a) shows an example of the identi- 



use in photonic applications, monolayer M0S2 transistors 
with a 10 8 on/off ratio and a room-temperature mobil- 
ity of more than 200 cm 2 /Vs have been demonstrated. 12 
This shows great potential in electronic applications ; 1 ^ 14 
since further improvements can be expected by, e.g. re- 
moving the ripples observed on monolayers* 1 ^ Although 
samples used in these studies&^ii^ were prepared by the 
simple mechanical exfoliation method^ techniques with 
potential for large-scale production, such as liquid and 
chemical exfoliations ^r— have been demonstrated to 
produce samples with high-mobility^ and similar PL 
properties^ 

For its applications in photonics and electron- 
ics, it is necessary to understand the dynamics of 
charge carriers in atomically-thin M0S2. Previ- 
ous steady-state optical studies, including PL ) 8 ' 9 ' 20 i 21 
absorption, 8,20 reflection, 8,9 photoconductivity, 8 and Ra- 
man scattering, 22 have revealed many aspects of elec- 
tronic and lattice properties. However, time-resolved op- 
tical measurements can provide direct information about 
the carrier dynamics, as illustrated by a recent time- 
resolved PL measurement^ 

Here we use a spatially and temporally resolved pump- 



fled flakes. From atomic-force microscopy measurements, 
one region has an average thickness of 1.5 nm. Since 
monolayer M0S2 (one Mo layer sandwiched by two S lay- 
ers) is 0.65-nm thick, we assign this region to a bilayer, 
considering uncertainties in the measurement. The re- 
gion next to it (to the left) has an average thickness of 
2.2 nm. It is assigned as a trilayer, which is next to a 
thick region of many layers (yellowish area) . A PL spec- 
trum obtained by focusing a 633-nm beam to the bilayer 
region is shown in Fig. 1(b). It peaks at about 675 nm. 
The sudden drop at the short-wavelength side is caused 
by a long-pass filter used to block the laser beam. Recent 
optical measurements have shown that the energy of the 
"A" exciton transition in atomically-thin samples is in 
the range of 655 - 685 nin tf 20 ' 23 Our result is consis- 
tent with these studies. When the excitation laser spot 
was moved to the trilayer region, no significant change in 
PL intensity was observed. This is a further confirmation 
that the first region is a bilayer instead of a monolayer, 
since PL from monolayer is expected to be much stronger 
than several-layer samples due to its direct bandgap. 8,9 

Figure [ljc) shows schematically the optical pump- 
probe setup. The pump pulse is focused to the sam- 
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FIG. 1. (a) Optical microscopy photo of the studied M0S2 
flakes under white light illumination, (b) PL spectrum of 
the sample, (c) Schematics of the optical pump-probe setup. 

(d) Schematics of the bandstructure of bulk M0S2 and the 
pump/probe configuration. The blue and red vertical arrows 
indicate the pump and probe photon energies, respectively. 

(e) Geometry of the multiple reflection from the sample. The 
rays are vertically shifted for clarity. 



pie surface by using a microscope objective lens. The 
probe pulse with a different wavelength is focused by the 
sample lens, and its reflection is collected and detected 
by a photodiode, the output of which is detected by a 
lock-in amplifier. The pump and the probe pulses are 
linearly polarized along perpendicular directions. A me- 
chanical chopper modulates the pump at about 2 kHz 
for lock-in detection. Reflection and scattering of the 
pump are prevented from reaching the photodiode by us- 

I 



ing color filters. In this configuration, the lock-in am- 
plifier measures the change in the reflection of the probe 
pulse from the sample induced by the pump. Such a 
quantity is then normalized to provide a differential re- 
flection, AR/Ro = (R— Rq)/Ro, where R and i?o are the 
reflection of the probe pulse from the sample with and 
without the presence of the pump pulse, respectively. All 
the measurements are carried out under ambient condi- 
tions and at room temperature. 

Figure 1(d) shows schematically the bandstructure of 
bulk M0S2. 8 Recent studies have shown that the indirect 
bandgap E' increases from the bulk value of 1.29 eV to 
more than 1.90 eV in monolayer, while the direct exci- 
tonic transitions (from VI and V2, respectively, to CI) at 
the K point remain largely unchanged. We use a pump 
pulse of 390 nm (3.179 eV) to excite the electrons from 
the valance band to the conduction band [the blue verti- 
cal arrow]. The pulse is obtained by frequency doubling 
a 100-fs, 780-nm pulse from a Thsapphire laser with a 
repetition rate of 81 MHz. The 660- nm probe pulse is 
obtained by frequency doubling the 100-fs, 1320-nm out- 
put of an optical parametric oscillator that is pumped 
by the Ti:sapphire laser. The probe photon energy of 
1.878 eV [the red vertical arrow] is near the "A" exciton 
transition from VI to CI at the K point. According to 
the PL spectrum shown in Fig. 1(b), the 660- nm probe 
is within the "A" exciton transition but higher than the 
peak energy. It is lower than the "B" exciton transition 
from V2 to CI of 2.04 eV (608 nm)l Hence, it detects 
the change of the " A" exciton transition induced by the 
carriers injected by the pump pulse. 

Figure 2(a) shows the measured differential reflection 
signal as a function of the probe delay, with a pump 
fluence of 11 /iJ/cm 2 . We find that the signal reaches 
a peak instantaneously (pulse-width limited) and exists 
for several hundred ps. The signal is positive, and is on 
the order of 10 _3 . The decay of the signal can be fit by 
a single exponential function [the red line in Fig. [U(a)], 
with a decay time of 100±10 ps. 

In order to use the differential reflection to monitor 
the charge carriers, it is necessary to establish the rela- 
tion between them. As shown in Fig. 1(e), the reflection 
beam is a superposition of three beams from the front 
surface and from the two interfaces. Hence, under nor- 
mal incidence, the reflection coefficient of this multilayer 
system can be written as22r— 
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where r x = (ho - hi)/(h + hi),r 2 = (hi - hi)/(hi + 
n 2 ),r 3 = (n 2 - h 3 )/(h 2 + h 3 ), and 4> t = 2Trd l n i /\. In 
these formula, hi, h 2l and h 3 are the complex indices 
of refraction for M0S2, Si02, and Si, respectively. The 



di are the thickness of the corresponding layers, and A is 
the wavelength. In general, the pump-injected carriers in 
M0S2 will change both the real part and the imaginary 
part of the index of refraction, h\ = n + ai, where n 
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FIG. 2. (a) Differential reflection signal measured with a 390- 
nm pump and a 660-nm probe. The solid line is a single 
exponential fit with a time constant of 100±10 ps. (b) Calcu- 
lated AR/Ro as a function of — Aa/ao by using Eq. 1, over a 
rather large range of AR/Ro- For small AR/Ro on the order 
of 10 -3 , the relation is approximately linear, (c) Calculated 
AR/Ro as a function of —An/no by using Eq. 1. 



and a are the real index of refraction and the absorption 
coefficient of M0S2, respectively. These will change the 
reflection coefficient according to Eq. I. 

Although the reflection coefficient depends on n and 
a in a rather complex way, for small changes, the dif- 
ferential reflection is expected to be proportional to the 
changes in n and a. We use Eq. I to verify this. By using 
published values of the n and a of thin- film MoS2^ we 
obtain i?o- We then vary a by a small quantity, Aa, and 
use Eq. 1 to calculate the corresponding AR. With this 
procedure, we obtain AR/Ro as a function of Aa/ao- 
The result is plotted in Fig. 2(b). Although a slight de- 
viation from a linear relation can be seen with relatively 
large Aa/ao, for small variations, AR/Ro is almost pro- 
portional to Aa/ao- Since the measured AR/Ro is on 
the order of 10 -3 , we can safely assume this linear rela- 
tion. By using a similar procedure, we verify the almost- 
linear relation between AR/Ro and An/no, as shown in 
Fig. 2(c). 

In our experiments, the probe pulse is tuned to the 
strong excitonic absorption peaks. Hence, we expect the 
absorption change to be the major contribution to the 
signal. That is, although Figs. 2(b) and 2(c) show that 
the same relative changes in n and a induce similar mag- 
nitudes of AR/Ro, the pump- injected carriers cause a 
larger Aa/ao than An/n - Because AR/Ro and Aa/ao 
have opposite signs, the positive signal shown in Fig. 2(a) 



indicates that the absorption of M0S2 is decreased by 
the pump-injected carriers. This is consistent with the 
phase space filling effect of free carriers and excitons^ 
which is often the dominant nonlinear effect observed in 
other materials like GaAs, although we could not exclude 
possible contributions from screening of the electron- 
hole interactional and bandgap renormalizationj 2 ^ For 
saturated absorption, a = a /(l — N/N sat ), where N 
and N sat are the carrier density and saturation density, 
respectively* 2 ^ Hence, Aa/ao — N/(N — N sat ). Since the 
observed Aa/ao is on the order of 10~ 3 , in these mea- 
surements N -C N sat and Aa/ao ^ N, approximately. 

With the established linear relation between the mea- 
sured AR/Ro and the carrier density, we can directly 
monitor the carrier dynamics. We can now attribute 
the 100-ps decay time obtained from the fit shown in 
Fig. [2fa) to the lifetime of the carriers excited by the 
pump pulse. We note that such a lifetime is consistent 
very well with the recently measured decay times of ex- 
citonic PL of about 100 ps^ 3 - 

The relatively long lifetime makes it possible to study 
diffusion of the carriers out of the pump laser spot dur- 
ing their lifetime. For a classical diffusion process with 
a gaussian initial profile that is determined by the pump 
laser spot, the squared width (full width at half max- 
ima) at time t is w (t) = w 2 (to) + 161n2£>(t - t ), where 
w 2 (to) is the squared width at an earlier time, and D is 
the diffusion coefficient. — Hence, by measuring the width 
of the carrier profile as a function of probe delay, we can 
deduce D. In this study, the produced flakes are rather 
small; hence, it is difficult to measure the full spatial pro- 
files. Therefore, we position the laser spots close to the 
boundary of a flake, and acquire one side of the profile by 
scanning the probe spot with respect to the pump spot. 
A few examples with several probe delays are shown in 
the inset of Fig. [3](a). We fit the measured profiles by 
gaussian functions, and the deduced widths are plotted 
in Fig. [31(b). Despite the large uncertainties, a trend of 
increase of the width with time is observed. By a linear 
fit [the red line in Fig. [3(b)], we deduce a diffusion coeffi- 
cient of 20±10 cm 2 /s. With the carrier lifetime of r = 100 
ps, this gives a diffusion length of Lp = v ' Dt = 450 nm. 

Owing to the Coulomb attraction, the excited elec- 
trons and holes move together in this diffusion process. 
Hence, the measured quantity is the ambipolar diffu- 
sion coefficient, D a — 2D e Dh/ (D e + Dh), where D e ^h) 
is the unipolar diffusion coefficient of electrons (holes) 4i 
Based on the fact that the effective masses of electrons 
and holes in M0S2 are similar^ 3 - we assume D e w Dh- 
Hence, the measured value can be treated approximately 
as unipolar diffusion coefficients of electrons and holes. 
Using the Einstein relation, D/kpT — /i/e, where fcs, 
T, /x, and e are the Boltzmann constant, the temper- 
ature, the mobility, and the absolute value of the elec- 
tron charge, the measured D corresponds to a mobility 
of fi w 800 cm 2 /Vs. Here, we have assumed the carriers 
have a thermal distribution with a temperature of 300 
K during the diffusion process. Since the measurement 
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FIG. 3. (a) Spatial profiles of the AR/Ro signal measured by 
scanning the probe spot with respect to the pump spot, with 
several probe delays as labeled in each panel. The solid lines 
are gaussian fits, (b) Square width of the spatial profiles de- 
duced from gaussian fits as a function of the probe delay. The 
solid line is a linear fit, corresponding to a diffusion coefficient 
of 20±10 cm 2 /s. 

was performed over a probe delay range of about 100 ps, 
which is much longer than typical thermalization and en- 
ergy relaxation time of a few ps in semiconductors, this 
assumption is valid. 

It is interesting to compare the optically measured mo- 
bility with electric measurements. Initial electric mea- 
surements on exfoliated monolayer M0S2 have indicated 
very low mobilities on the order of 1 cm 2 /Vsi£. Very re- 
cently, it was found that by covering the M0S2 layer with 
a dielectric layer with a high dielectric constant, such as 
Hf02, the mobility can be increased to 780 cm 2 /Vs for 
monolayer MoS2- 12 This was attributed to suppression of 
Coulomb scattering by the dielectri o 12 ' 32 ' 33 and modifica- 
tion of phonon dispersion ! 12 ' 34 The M0S2 flakes studied 
here are not covered by a dielectric layer. However, the 



mobility deduced is comparable to those with a dielectric 
top layer. We attribute this fact to the ambipolar nature 
of the diffusion process: In this process, the electron- 
hole pairs move as a unit. Since the pair is electrically 
neutral, it is less influenced by the Coulomb scattering. 
Therefore, although the M0S2 sample is not covered by 
a dielectric layer, the mobility we measured indicates to 
the one that is not limited by the Coulomb scattering. 
We note that such a comparison is also limited by the 
exciton formation from the electron-hole pairs, which is 
enhanced in atomically-thin layers due to the quantum 
confinement. 

In summary, we have used a femtosecond pump-probe 
technique to study charge carriers in few-layer M0S2 sam- 
ples. Carriers are injected by direct interband absorp- 
tion of a 390-nm pump pulse, and detected by measur- 
ing the differential reflection of a probe pulse of 660-nm. 
We found that the absorption is reduced by the carriers, 
which is consistent with the phase-state filling effect. By 
spatially and temporally resolving the differential reflec- 
tion signal, we deduce a carrier lifetime of 100±10 ps and 
a carrier diffusion coefficient of 20±10 cm 2 /s, correspond- 
ing to a mobility of 800 cm 2 /Vs and a diffusion length of 
450 nm. Besides providing quantitative information on 
the carrier dynamics in this new promising nanomaterial, 
our experiment may stimulate further optical studies of 
carrier dynamics in this material system. The optical ap- 
proach has the advantage of being noninvasive, avoiding 
device fabrication and influence of electric contacts. 
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